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Abstract 

A linear nonadiabatic analysis of the vibrational stability of population III main-sequence stars was 
carried out. It was demonstrated that, in the case of massive stars with M ^ 5Mq, helium burning (triple 
alpha reaction) starts during the main-sequence stage and produces ^^C, leading to the activation of a 
part of the CNO-cycle. It was found that, despite of that, those stars with M ^ 13Mq become unstable 
against the dipole gi- and g2-modcs due to the e-mcchanism, during the early evolutionary phase at which 
the pp-chain is still the dominant nuclear energy source. The instability due to the e-mechanism occurs 
against g-modes having a large amplitude in the off-centered '^He accumulation shell in the deep interior, 
and the growth time is much shorter than the evolutionary timescale. This instability is therefore likely to 
induce mixing in the stellar interior, and have a significant influence on the evolution of these stars. 
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1. Introduction 

The structure and evolution of the first-generation stars 
(so-called population III stars) born in the very early uni- 
verse with no, or little, heavy elements are considerably 
different from those of stars born much later with heavy 
elements (e.g., Weiss et al. 2000, Marigo et al. 2001, Suda 
et al. 2007). This is partly because the initial mass func- 
tion is substantially different due to the lack of an efficient 
coolant at the star-formation stage, and also partly be- 
cause the only possible way of starting of hydrogen burn- 
ing in population III stars is the pp-chain, of which tem- 
perature dependence is weak, no matter how massive the 
star is. Hence, the temperature near the stellar center 
of intermediate and massive population III stars is much 
higher than in the case of the population I/II stars, and 
then these population III stars are, to a great extent, com- 
pact compared with the population I/II stars. Since the 
pp-chain is the main energy source, the temperature gra- 
dient in the nuclear burning core of mctal-frcc population 
III stars is milder than in the case of population I stars 
with heavy elements. Consequently, the nuclear burning 
core is radiative, except for during the short period of the 
ZAMS stage. Furthermore, the high central temperature 
of population III stars leads to the occurrence of helium 
nuclear burning (triple alpha reaction) even during the 
main-sequence stage. On the other hand, no metallicity 
and the consequent low opacity make the envelope to be 
radiative, except for the surface layer. 

Since the structure of population III stars is consider- 
ably different in this way from that of stars with heavy 
elements, their stability nature might also be different. 



The present authors have already shown by using a quasi- 
adiabatic approximation that low-mass population III 
stars with M ^ 5 Mq become unstable during the main- 
sequence stage against the dipole low-order g-modes due 
to the e-mechanism (Sonoi & Shibahashi 2011, hereafter 
Paper I) . The physical cause of this instability is explained 
as follows (Christenscn-Dalsgaard et al. 1974; Boury et al. 
1975; Shibahashi et al. 1975; see also Paper I): ^He accu- 
mulates in a shell just outside the stellar center in the 
process of the pp-chain, since the '^He-'^He reaction occurs 
so fast compared with the pp-reaction that '^He does not 
remain at the stellar center, while the ^He-^He reaction 
does not efficiently occur in the outer part of the nuclear 
burning core with a decrease in temperature. In such a 
situation, the high-temperature sensitivity of the •^He-^He 
reaction leads to a vibrational instability of g-modes hav- 
ing a large amplitude of temperature perturbation in such 
an off-centered '^He shell. 

In Paper I, when examining the vibrational stability of 
the g-modes, we adopted the quasi-adiabatic approxima- 
tion, in which the work integral was evaluated by using 
the adiabatic eigenfunctions as well as adiabatic relations. 
This approximation is good enough in the bulk of the stel- 
lar interior, but ibecomes inappropriate in the very outer 
envelope where the thermal time scale is short and the 
nonadiabatic effect is large. To avoid this difficulty, the 
integration was terminated at a radius where the ther- 
mal timescale became equal to the oscillation period. As 
a consequence, the contribution of the K-mechanism was 
not properly evaluated, though there appeared a tendency 
for a destabilization effect due to the K-mechanism in some 
modes. In the present work, we carried out a fully non- 
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adaiabatic analysis to see the additional effect of the k- 
mechanism on the vibrational stability of the g-modes. 

In Paper I, we restricted ourselves to treat only low- 
mass stars, M < 5Mq, in which the CNO-cycle never 
starts to work at the early stage of core hydrogen-burning. 
In the case of more massive stars, even at the early stage of 
core hydrogen-burning, the central temperature becomes 
high enough for helium burning (triple alpha reaction) 
to start in addition to hydrogen burning due to the pp- 
chain. In such a case, the production of ^^C by the triple 
alpha reaction leads to activation of, at least, a part of 
the CNO-cyclc, and hence a careful treatment of the nu- 
clear reaction networks is required. In the present work, 
we carried out stellar evolution calculations with nuclear 
networks, and extended the mass range of stars. 

The purpose of the present work is to examine the vi- 
brational stability of population III stars more completely 
than Paper I by performing stability analysis by solving 
the fully nonadiabatic oscillations over a wide mass range 
of population III stars. 

2. Evolutionary equilibrium models of popula- 
tion III stars 

We adopted a recently developed, open source stellar 
evolution code, MESA (revision 3107, Paxton ct al. 2011), 
which is appropriate for the current purpose of provid- 
ing state-of-the-art modules for equation of state, opacity, 
nuclear reaction rates, element diffusion data, and atmo- 
sphere boundary conditions. Evolutionary models of 1, 
1.2, 1.4, 1.6, 2, 3, 5, 7, 10, 13, 16, and 2OM0 with the ini- 
tial chemical composition Xq = 0.76 and Z — were con- 
structed from the pre-main sequence Hayashi phase. Since 
this code carefully treats the change in the abundance of 
chemical elements, the evolution of population III stars is 
computed with sufficient precision even after the onset of 
the triple alpha reaction and the resultant activation of 
the CNO-cycle. On the other hand, this code outputs, at 
default, physical quantities only below the photosphere; 
we therefore modified it to output an Eddington gray at- 
mosphere as well by utilizing the MESA's module so that 
we could set the outer boundary conditions of the pulsa- 
tion equations at an optically thin layer (From revision 
3372, MESA has given an option to output atmosphere 
models.). 

A lack of heavy elements leads to low opacity, and 
thus the luminosity of population III stars is substantially 
higher than in the case of stars with heavy elements. This 
means that the population III stars are compact com- 
pared with population I stars; also, the central density 
and temperature are significantly higher than in the case 
of population I stars. This feature is clearly seen in fig- 
ure 1. A conspicuous difference between metal- free stars 
and stars with Z ^ in this diagram is that the central 
density becomes higher and higher with an increase in 
the mass of the star, up to M ~ I6M0, in the case of 
metal-free stars, while in the case of Z 7^ it turns to be- 
come lower with an increase of the mass of the stars, for 
M > 1.2 Mq. This difference is caused by the fact that 
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Fig. 1. Physical conditions at the center of 1, 1.2, 1.4, 1.6, 
2, 3, 5, 7, 10, 13, 16, 20, 4OM0 stars at the ZAMS stage with 
different metaUicity. The abscissa is the density in units of 
kgm""^, and the ordinate is the temperature. 

the energy source of population I/II stars turns to be the 
CNO-cycle for M > 1.2 Af© from the pp-chain, while in the 
case of metal-free stars the pp-chain remains as the only 
source, even in the mass range of M ^ 1.2M0. While there 
appears to be a convective core in the case of metal-free 
stars with I.2M0 < M < 16Mq at the ZAMS stage, and 
very close to it, despite that the nuclear reaction source is 
the pp-chain, due to the high-temperature gradient, the 
nuclear burning core soon becomes radiative in the case of 
lower mass. On the other hand, for stars with M > 16 M©, 
the central temperature becomes high enough for helium 
burning to occur at the ZAMS stage. As a consequence, 
the CNO-cycle is activated with the produced ^^C at zero- 
age, and this process dominates over the pp-chain. This 
is the reason why the central density becomes lower with 
an increase of the mass of the star for M ^ 16 M©. The 
nuclear burning core of these stars is convective. Figure 2 
shows the evolutionary tracks of the stars on the HR dia- 
gram. Since the nuclear reaction occurring in population 
III stars of M< I6M0 at the ZAMS stage is the pp-chain, 
they go toward higher temperature and higher luminosity 
on the HR diagram, like in the case of M < 1.2 Af© of 
population I stars. 

As can be seen in figure 3, in the case of 2M0 ^ Af ^ 
16 A/0, though the contribution of helium burning to the 
total luminosity is negligibly smaller compared to that of 
the pp-chain, it gradually becomes large enough during 
the main-sequence stage to produce ^^C, which activates 
the CNO-cycle of hydrogen burning (Suda et al. 2007). 
The CNO-cycle onset stage is marked by open squares in 
figure 2. Since the temperature dependence of the CNO- 
cycle is much higher than that of the pp-chain, the CNO- 
cycle becomes to dominate over the pp-chain if once the 
CNO-cycle is activated. Also, the nuclear burning core be- 
comes convective. Figure 3 shows the evolutionary change 
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Fig. 2. Theoretical HR diagram showing the evolution- 
ary tracks of 1.0, 1.2, 1.4, 1.6, 2.0, 3.0, 5.0, 7.0, 10.0, 
13.0, 16.0, and 20.0 Mq stars with the initial composition 
(X, Z) = (0.76, 0), (0.70, 0.02). The bold lines on the Z = 
tracks indicate stages from Xc = 0.76 to 0.01 {Xc is the cen- 
tral hydrogen abundunce), in which the stability analysis is 
performed in this study. In the case of M> 1.27^0 stars, evo- 
lutionary tracks of metal-free stars are significantly different 
from those o{ Z j^O because of suppression of the CNO-cycle. 

in the contribution of the CNO-cycle to the total nuclear 
energy generation for stars wherer M > 2 Mq . As can be 
seen in this figure, the contribution of the CNO-cycle in- 
duced by triple alpha reaction becomes more important 
with an increase of the mass of the star. For stars with 
M ^ 16 Mq, the CNO-cycle works even at the zero-age. 
On the other hand, for stars with M ^ 2Mq, the nuclear 
energy is provided only by the pp-chain in most of core 
hydrogen-burning phase. 

The characteristic properties of a population III star as 
an oscillator have already been discussed in Paper I. It is 
remarked here again that, for a given mass, the radius of 
the population III star is substantially less than in the case 
of population I stars, and hence the dynamical timescale, 
which gives a rough estimate of the oscillation periods, 
of population III stars is substantially shorter than that 
of population I stars. Another important factor for the 
oscillation periods of g-modes is the profile of the Brunt- 
Vaisala frequency. In the case of low-mass population III 
stars, an absence of the convective envelope increases the 
g-mode frequencies in units of the dynamical timescale 
of the star. Combining these two factors, the oscillation 
periods of the g-modcs in population III stars are shorter 
than in the case of population I stars, as long as the mass 
and the mode are fixed. 

3. Fully nonadiabatic analysis of vibrational sta- 
bility 

Any eigenfunction of linear nonradial oscillations of 
stars is expressed in terms of the combination of a spatial 
function and a time-varying function. The latter is ex- 
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Fig. 3. Fraction of nuclear energy by the CNO-cyclcs and 
triple alpha reaction to the total nuclear energy at the stellar 
center. The abscissa is the central hydrogen abundance, Xc, 
which is an indicator of stellar evolution. This range is set 
to be from X^ = 0.76 to 0.01. For < 1.6 Mq, the nuclear 
energy is provided almost exclusively by the pp-chain. The 
contribution of the triple alpha reaction to the whole nuclear 
energy is negligible. 

pressed by exp(icrt), where a denotes the eigenfrcquency. 
The former, the spatial part, is decomposed into a spheri- 
cal harmonic function, which is a function of the colatitude 
and the azimuthal angle, and a radial function. Equations 
governing the radial functions lead to a set of sixth-order 
differential equations, of which coefficients are complex, 
including the terms of frequency a. Together with the 
proper boundary conditions, this set of equations forms 
a complex eigenvalue problem with an eigenvalue a. The 
real part of a, ctr,, represents the oscillation frequency, 
and the imaginary part of tr, cti, gives the growth rate, 
or the damping rate, depending on its sign. We have 
newly developed a Henyey-type code to solve this eigen- 
value problem, following the formulation by Unno et al. 
(1989). 

For a stability analysis relevant to the e-mechanism, it 
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is important to evaluate the temperature and density de- 
pendence of nuclear energy generation on the oscillation 
time scale, et = {d\ne/d\nT)p and Ep = (91ne/51np)T- 
It should be noted here that wc should adopt the depen- 
dences of the nuclear reaction through perturbations, and 
that they are different from those for the evolutionary 
time scale. For example, as for the pp-chain, the tem- 
perature dependence of energy generation all through the 
pp-chain in equilibrium is governed by the slowest reac- 
tion, ^H(^H,e+i'o)^H, and et — 4; we should thus adopt 
the effective temperature dependence of the nuclear reac- 
tion through the perturbation, which is mainly governed 
by 3He(3He,2iH)4He and et ^ H (Dilke & Gough 1972; 
Boury & Noels 1973; Unno 1975; Unno et al. 1989). We 
evaluated these effective dependences of the pp-chain and 
of the CNO-cycle through the perturbation, separately, 
and then averaged them with their contribution to the 
total nuclear energy generation to obtain the net values. 

The temperature and density dependence of the pp- 
chain were evaluated following Unno et al. (1989). For 
the CNO-cycle we took only the CN-cycle into account, 
since the occurrence probability of the other processes 
of the CNO-cycle is very low (~ 10"*). The CN-cycle 
involves relatively slow /3-decay, of which the lifetime 
could be comparable to the oscillation period. Kawaler 
(1988) derived the formulation for the temperature and 
density dependence, while taking into account the phase 
lag between the creation and destruction of various re- 
actants. He used an approximation assuming that the 
time scales of /3-decays in the CN-cycle ^^N(e~,i/e)^'^C and 
^^0(e~, i'e)^^N are comparable to the oscillation period 
n, while those of the other collisional reactions are much 
longer than IT. However, since in population HI stars the 
onset of the CNO-cycle occurs at extremely high temper- 
ature lO^K), just after ^^C is produced by the triple 
a reaction, the time scale of the collisional reactions can 
also be comparable to the oscillation period H. Therefore, 
we evaluated the temperature and density dependences 
of CN-cycle without assuming that the collisional reac- 
tions are much slower than the oscillation. For the sake of 
simplicity, we ig nored the branch ^3N(p,7)"0(e~,i/e)"N, 
which must occur at ^ 10^ K. The details of the evalua- 
tion of the density and temperature dependence of the 
CN-cycle are given in the Appendix. 

The stability is directly determined by the imaginary 
part of a. It should be noted here that the work integral 
is useful to see where the excitation and damping zones 
are located in the stellar interior. The work integral, W, 
is defined as the amount of energy that must be removed 
from the star in order for the star to oscillate strictly pe- 
riodically, and it is related to the imaginary part of cti 
through 
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where H — 27rcrj^ denotes the period and Ew is the total 
energy of oscillation, which is twice the time average of 
the oscillation kinetic energy. 
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where ^ denotes the eigenfunction of the displacement. 
The work integral was evaluated in Paper I using the 
quasi-adiabatic approximation. In the present study, we 
carried out fully nonadiabatic calculations, and evaluated 
the work integral by 
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where 5T and 5S denote the eigenfunction of the 
Lagrangian perturbation of temperature and that of en- 
tropy, respectively. W{r) represents the work done on the 
overlying layer by a sphere of radius r; its surface value is 
W, used in equation (1). 

4. Instability due to e-mechanism of the pp-chain 

As discussed in section 1, the favorable modes for the 
£-mechanism of '^He-'^He reaction are the low-degree g- 
modes having a large amplitude in the off-centered '^He 
accumulation shell, which is near to the outer edge of the 
nuclear burning core. Therefore, we investigate the dipole 
(/ = 1) and quadrupole {I = 2) modes. 
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Fig. 4. Work integrals for diploic {I = 1) gi-mode in stars of 
1-20 Mq . Equilibrium models correspond to the filled circles 
in the top panel of figure 5. 

Figure 4 shows the work integral for the dipole gi- 
modes for the models of 1, 1.2, 1.4, 1.6, 2, 3, 5, 7, 10, 
13, 16, and 20 M© stars with the initial hydrogen abun- 
dance Xq = 0.76 and Z = 0. Each model corresponds to the 
most unstable (or the least stable) model along the evo- 
lutionary stage at around Xc — 0.6, which is marked as a 
filled circle in the top panel of figure 5. As can be clearly 
seen in this figure, the nuclear burning core contributes to 
dcstabilization through the e-mechanism, while the bulk 
of the outer part, except for the zone around logT ~ 5.2 
in the case of M ^ 2Mq, works against it. The positive 
contribution from logT ~ 5.2 is the K-mechanism of the 
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the second helium ionization zone, which is discussed later 
in detail. 

As can be seen in figure 3, the main nuclear source of 
stars with M ^ 13 Mq around at the evolutionary phase 
of Xc — 0.6, which are shown to be vibrationaly unstable, 
is the pp-chain. On the other hand, in the case of more 
massive stars, M ^ 16 Mq, which are shown to be vibra- 
tionaly stable, the main nuclear source is the CNO-cycle. 
As a consequence, we reach the conclusion that the popu- 
lation III stars with M ^ 13M0 are unstable against the 
dipolc gi-modc, due to the e-mcchanism of the pp-chain. 
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Fig. 5. Top: Variation of tlie growtii rates of the dipole 
{I = 1) gi-mode witii stellar evolution. The filled circles corre- 
spond to the equilibrium models of which the work integrals 
arc shown in figure 4. Bottom: Same as the top panel, but 
for the dipolc g2-modc. 

Figure 5 shows the growth rates of the dipole gi- and 
g2-modes as a function of the hydrogen mass function at 
the center, Xc, which is an indicator of stellar evolution. 
The dipole gi-mode becomes unstable from Xc — 0.7 to 
0.4 for 1.2Mq <M < IOMq. For IM©, instability grows 
rapidly after Xc — 0.4. We discuss it later in the next 
section. 



The growth rates of the gi- and g2-modes increase due 
to the £-mcchanism of the pp-chain with stellar evolution 
after leaving the ZAMS (figure 5). As explained in Paper 
I, this delicate change of stability is caused by the varia- 
tion in eigenfunction with evolution. A high-temperature 
sensitivity of the nuclear energy generation rate of the 
'^He-'^He reaction is essential for the £-mechanism to work 
in the case of the lower mass stars (Christensen-Dalsgaard 
et al. 1974). The effective value of ex in the perturba- 
tion to the nuclear energy generation is raised to 11 
because of •^He("^He, 2^H)^He, while the value of Et of the 
pp-reaction, itself, ^H(^H,e"''i'o)^H; is ~ 4 ~ 5. Since the 
reaction '^He('^He, 2^H)^IIe occurs fast compared with the 
pp-reaction, '^He does not remain at the stellar center. On 
the other hand, the '^He-'^He reaction does not efficiently 
occur with a decrease in temperature. As a consequence, 
■^He accumulates in a shell just outside the stellar center. 
Hence the most favorable situation for the vibrational in- 
stability is that the eigenfunction has a large amplitude 
in such an off-centered '^He shell. At the ZAMS stage, the 
accumulation of ^He is not enough, and thus the star is 
stable, or marginal, against any modes. As the evolution 
proceeds, the peak of the eigenfunction becomes to match 
with the off-centered '^He shell, and then the star becomes 
vibrationally unstable against the low-order g-modes. For 
the higher mass stars, on the other hand, the pp-II and III 
branches contribute to the nuclear energy generation due 
to high temperature. The collisional reactions in these 
branches have high temperature dependence comparable 
with that of the '^He-'^He reaction, and contribute to the 
g-mode instability. 

In the case of stars with M > 13 Mq, the CNO-cycle is 
activated to work as the dominant source of the nuclear 
reaction at the evolutionary stage, Xc — 0.7. This is the 
reason why the unstable duration for the star with M = 
13AfQ is much shorter than for the less-massive stars, and 
why the more-massive stars do not have instability due to 
the e- mechanism (see the top panel of figure 5). 

In the case of a 10 star, after the phase of Xc — 0.4, 
the growth rate changes rapidly, and it grows to be posi- 
tive in the case of the g2-mode. This is because the CNO- 
cycle becomes the dominant source of energy generation, 
and as a consequence the convective core grows without 
vanishing. A superadiabatic zone then appears outside 
the convective core, and some g-modes become excited by 
Kato's mechanism (Kato 1966). We will discuss this insta- 
bility and the relevant physics in detail in the forthcoming 
paper (in preparation). 

5. The K-mechanism as a booster for instability 

It is well known that the K-mechanism in the hydrogen 
and helium ionization zones is responsible for the Cepheid 
strip of radial pulsations. The depth of the helium ioniza- 
tion zone in stars having logToff = 3.8 ^ 4.0 is appropriate 
to excite pulsation modes, of which frequencies arc on the 
order of {GAI/R^y^^ by the K-mechanism. The S Set 
stars, which arc around M ~ 2Mq, are such pulsating 
stars near to the main sequence in the case of population 
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I stars. 

Since the opacity is much lower in the metal-free popu- 
lation III stars than in the case of population I stars, the 
luminosity of population III stars is higher than that of 
population I stars of the same mass. Also, the population 
III stars are more compact than the population I stars. As 
a consequence, the effective temperature of population III 
stars is substantially higher than the population I stars. 
As cam be seen in figure 2, in the case of metal- free pop- 
ulation III stars, the main-sequence stars of M ~ 1 Mq is 
around log Toff = 3.8 ^ 4, and hence the K-mcchanism of 
the the second helium ionization zone does work to excite 
low-degree, low-order p-modes in these stars. 
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Fig. 6. (a): Work integral for the dipole pi mode of a 
1 Mq star model at the evolutionary stage of Xc = 0.71, 
logTefi = 3.86. The solid line shows the results obtained by 
the fully nonadiabatic calculation, while the dashed line shows 
that obtained by the quasi-adiabatic approximation. The cir- 
cle on the dashed curve for the quasi-adiabatic approximation 
indicates the location where the oscillation period is equal to 
the local thermal time scale, (b): Same as (a), but for the 
dipole gi- and g2-modes. (c): Same as (a), but for the dipole 
gi2 mode, (d): Opacity, k, the nuclear reaction rate, e, and 
the fraction of the radiative luminosity to the total luminosity, 
Lp/L, as functions of logT. 

Panel (a) of figure 6 shows the work integral for the 
dipole pi-mode. As can be seen in this figure, the p- 
modes are excited by the K-mechanism, mainly due to 
the first helium ionization. Many of the other higher or- 
der p-modes and higher degree p-modes are also unstable. 
This situation is analogous to the situation of population I 
stars in the classical instability strip (Shibahashi & Osaki 
1981). It should be noted that, since the p-modes os- 
cillate mainly in the outer envelope, where the density 
is much lower than in the deep interior, their amplitude 
grows more easily than g-modes, which oscillate in the 
deep interior. Hence, the growth rate of the modes ex- 
cited by the K-mechanism is generally larger than those 
excited by the e-mechanism. 

The positive contribution from logT ~ 5.2 seen in fig- 
ure 4 in the case of M ^ 2 Mq is the K-mechanism of the 
the second helium ionization zone. Panel (b) of figure 6 
shows the work integrals of the dipole gi- and g2-modes 



for a IMq star model with log Toff = 3.86. This demon- 
strates that the K-mechanism boosts up the work integral 
at the second helium ionization zone at around logT~5.2, 
and the first helium ionization zone at around logT ~ 4.7 
after the contribution of e-mechanism is almost cancelled 
by the radiative dissipation outside the nuclear burning 
core. It should be reminded here that in the case of the 
quasi-adiabatic approximation we cannot properly eval- 
uate the contribution of the very outer envelope where 
logT ^ 4.7 to the work integral, and we have to terminate 
the integration around the zone where the local thermal 
timescale becomes as short as the oscillation period. The 
current fully nonadiabatic calculation shows that the heat 
capacity above such zones is so small that the contribu- 
tion of the first helium ionization zone and the hydrogen 
ionization zone are negligibly small. 

As the evolution proceeds, the Brunt- Vaisala frequency 
in the stellar interior increases due to an increase of the 
central condensation and the formation of the gradient of 
chemical compositions as a consequence of nuclear reac- 
tions converting hydrogen to helium, and hence the fre- 
quencies of g-modes increase. When the frequency of a 
g-mode is close to a p-mode, these two modes become to 
have a dual character: behaving with a g-mode character 
in the deep interior while with a p-mode character in the 
envelope. These eigenmodes never degenerate, however, 
and avoided crossings of modes occur, — that is, the wave 
characters are switched for each other (see Unno et al. 
1989). 

Figure 7 (a) shows the evolutionary variation in the 
nondimensional frequency, wr = aB./{GM/R^y/^ , nor- 
malized with the dynamical time scale, and the stability 
of the dipole and quadrupole modes of the 1 Mq star. In 
this diagram, the small dots denote the stable modes. As 
for the unstable modes, we classify them into four types 
by the degree of contribution of the £-mechanism to the 
work integral: (i) the modes where the e-mechanism con- 
tributes more than 50%, shown with filled squares, (ii) 
those with 30 — 50%, shown with filled circles, (iii) those 
with 10 — 30%, shown with filled triangles, and (iv) those 
with less than 10%, shown by open circles. The modes 
of the type (i) can be regarded as being excited mainly 
by the e-mechanism, while the type (iv) are regarded as 
being excited essentially by the K-mechanism. 

As can be seen in this figure, the p-modes of the 1 Mq 
star are vibrationally unstable during the main sequence 
stage. This is because the effective temperature of the 
metal-free 1 Mq star is high enough to be in the range of 
the instability strip on the HR diagram. Those p-modes 
are excited by the K-mechanism working in the second he- 
lium ionization zone. With stellar evolution, the frequency 
of the gi-mode approaches that of the pi-mode, and the 
gi-mode gradually becomes to have a mixed character of 
a p-mode and a g-mode. As a consequence, the relative 
contribution of the K-mcchanism to the work integral be- 
comes larger and larger than that of the e-mechanism. 
Hence, the growth rate of the gi-modc becomes signifi- 
cantly larger. This is the reason why the growth rate of 
the dipole gi-mode significantly increases. 
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Fig. 7. (a); Variation in the frequencies and stability nature 
of tile dipole {I = 1) (upper panel) and quadrupole {I = 2) 
(lower panel) modes for the metal-free 1 Mq star with evo- 
lution. The small dots denote stable modes, while the big- 
ger symbols denote unstable modes. The latter are catego- 
rized into four types: (i) the modes that contribution of the 
£-mechanism to the work integral is more than 50%, shown 
with the filled squares, (ii) those that contribution of the 
e-mechanism is 30-50%, shown with the filled circles, (iii) 
those that contribution of the e-mechanism is 10-30%, shown 
with the filled triangles, and (iv) those excited essentially by 
the K-mechanism, shown with the open circles. The abscissa 
is the central hydrogen abundance Xc, which is an indica- 
tor of stellar evolution, (b); Same as (a), but for the dipole 
modes for the 1.2 Mq star (upper panel) and for the 1.4Mq 
star (lower panel). 

With the increase of the stellar mass, the effeetive tem- 
perature at the main-sequence stage becomes higher. The 
K-mechanism then becomes less efhcient to excite oscil- 
lations. Marigo et al. (2001) showed that the instability 
strips for the radial fundamental mode, first and second 
overtones cross the main-sequence track for the popula- 
tion III ~ IMq star. Figure 7 (b) shows that in the case 
of a 1.2 Mq star, the dipole pi-mode, of which the char- 
acteristics are similar to the radial fundamental mode, is 
stable during the hydrogen-burning phase, and that also 



the other p-modes are stable, except at the very early 
phase of evolution. In the case of I.AMq, the p-modes are 
no longer excited. 

6. Stability of low-degree, high-order g-modes 

It is found the K-mechanism does work to excite low- 
degree high-order g-modes as well as p-modes simultane- 
ously during the early phase of evolution. Panel (c) of 
figure 6 shows the work integral for the dipole gi2-mode. 
The amplitude of such a high-order g mode in the '^He 
burning shell is relatively much lower than in the outer 
envelope, and thus the contribution of the £-mechanism 
to the stability of this mode is negligibly small. Rather, 
this mode is excited by the K-mechanism of hydrogen ion- 
ization. The stability is, however, marginal. 

One might be reminded of hybrid population I pulsators 
of M ~ 1.5 Mq, expected theoretically in relation to 6 Set 
and 7 Dor variables. Guzik et al. (2000) first proposed 
"convective blocking" as the driving mechanism of high- 
order g-modes in the 7 Dor variables. Dupret et al. (2004, 
2005) performed a stability analysis while adopting the 
time-dependent convection theory. They found that this 
mechanism is dominant for the excitation, and hence that 
perturbation of convective energy flux has no significant 
effect. However, the situation is substantially different. 
In the present case of population III stars, the convective 
layer is induced only by hydrogen ionization, and it is 
much thinner than the population I stars; also, the high- 
order g-modes are excited by the K-mechanism near to the 
top of the convective layer. Since the convective timescale 
there is shorter than the pulsation timescale, it is still hard 
to conclude definitely the stability of such case. 

7. Summary and discussion 

We performed a nonadiabatic analysis of the vibrational 
instability of population III stars due to the £-mechanism 
of the pp-chain, which was analyzed with a quasi-adiabatic 
approximation in Paper I for stars M < 5Mq. We ex- 
tended the mass range by adopting a stellar evolution 
code, MESA, which treats more appropriately the nu- 
clear reaction networks after the onset of the CNO-cycle. 
We found that the dipole gi-mode becomes unstable in a 
wider range of masses than investigated in Paper I, and 
that the stars of M ^ 13 Mq become unstable at their early 
evolutionary phase, in which the pp-chain is the dominant 
energy source. It was found that the dipole g2-mode also 
becomes vibrationally unstable for some low-mass stars. 
This mode, however, is found to be unstable for a quite 
short fraction of the star's lifetime. This mode becomes 
unstable for stars of M < I.2M0 and 3Mq <M < 5Mq. 

The e-folding time of the instability of the dipole gi- 
mode is much shorter than the evolutionary timescale. 
Hence, it seems to be natural to expect that the instability 
grows to a finite amplitude. If the oscillatory wave grows 
to a huge amplitude, it is expected to induce material 
mixing, and to have a significant influence on the later 
evolution of the star, as was expected concerning the Sun 



8 



T. Sonoi and H. Sliibahashi 



by Dilke & Gough (1972). 

One of the possible influences on the stellar evolution 
might be an extension of the lifetime. That is, incorpora- 
tion the surrounding cooler and hydrogen-rich matter into 
the central region could delay the hydrogen exhaustion. It 
has already been expected that the lifetimes of the popu- 
lation III stars with < O.SAf© should be longer than the 
age of the Universe, and should still be alive if ever born 
(Marigo et al. 2001). If material mixing could be induced 
by the instability, the mass limit for the population III 
stars being still alive would be increased. Hence, a pos- 
sibility of detecting low-mass population III stars might 
be higher than expected before, and there might remain 
a possibility of finding oscillations in such stars by moni- 
toring a micro-magnitude level luminosity variation. 

Population III stars with M ~ IM© were found to be 
unstable for both of the dipole g-modes and the low-degree 
p-modes simultaneously. The g-modes are excited by the 
e-mechanism of the ■^Hc-'^He reaction, and the p-modcs are 
excited by the K-mcchanism in the helium ionization zone. 
This is a quite unique hybrid pulsator, only expected in 
metal-free population III stars. In the case of population 
I stars, the main-sequence stars in the classical instability 
strip are ~ 2Mq stars, in which the nuclear energy source 
is the CNO cycle, and their g-modes arc not excited by the 
e-mechanism. At a certain evolutionary phase, avoided 
crossing occurs between the g-modes and p-modes. In 
such a case, the e-folding time of the g-modes significantly 
increases with the help of the K-mechanism and becomes 
substantially shorter than the evolution timescale. 

This study has concentrated on the main-sequence 
stars. It should be noted here that the number of known 
red giants with extremely mcta-poor abundances is larger 
than that of dwarfs by a factor of two (Suda et al. 2008). 
Also, it should be reminded that solar-like p-mode oscil- 
lations in red giants with normal abundances are being 
detected in tens of thousands of stars with space-based 
photometric campaigns. Hence, one might be able to ex- 
pect detection of solar-like p-mode oscillations in popula- 
tion HI giants in the future. It should be remarked here 
that the extent of convectivc envelope might be different 
between population I stars and population HI stars due to 
difference in the opacity. Thus, the convection spectrum, 
and hence property of stochastic excitation, are also ex- 
pected to be different. 

This study has been financed by Global COE Program 
"the Physical Sciences Frontier", MEXT, Japan. It 
adopted the MESA code to construct the stellar evolu- 
tionary models. The authors wish to acknowledge the 
Paxton et al. (2011) publication and MESA website^ for 
providing their very useful evolution code to public as an 
open source. 
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Appendix. Evaluation of the density and temper- 
ature dependence of nuclear energy generation by 
the CN-cycle 



Table 1. The labels for the elements in the CN-cycle. 



i 


species 


reaction 


1 




^^C(p,7)^-^N 


2 


13N 


i3N(e-,z/e)"C 


3 




"C(p,7)"N 


4 




i%(p,7)i50 


5 


150 


i50(e-,i/e)^^N 


6 


15N 


i5N(p,4He)i2C 



This section introduces how to evaluate the density and 
temperature dependence of the nuclear energy generation 
by the CN-cycle, £p,cN and £t,cn- We slightly modified 
the method shown in Kawaler (1988) to obtain a reason- 
able result applicable even in logT ^ 8. 

If we label the species of the CN-cycle and proton as 
i's, following table 1 and p, respectively, the differential 
equations governing the abundances of the rcactants in 
the CN-cycle are written as follows: 

dt 
dt 
dt 

where TV denotes the total number density and yi the num- 
ber fraction of the element i; Ci denotes the rate of coUi- 
sional destruction of species i per proton for i ~ 1,3,4,6; 
and Xi is the rate of decay of species i for z = 2,5. Here, the 
label J means i — I, except for the case of z = 1 , in which 
the label j on the right-hand-side of equation (Al) should 
be recognized as ^^N (z = 6). The Lagrangian variation of 
j/i's is expressed as 



N^y.VpC^ 4- 




(i = l,4) 


(Al) 


Ny,X + N^ 


j/i-ij/pCi-i, 


(z = 2,5) 


(A2) 


N^y^ypC^ + 


7V?/,_iA,_i, 


(z = 3,6) 


(A3) 



Syi Sp 6T 
— = l-ftT77- 



(A4) 



That is, the coefficients a^'s and /3i's are the density 
and temperature dependence of yi, respectively. By us- 
ing these coefficients, the temperature and density depen- 
dences of the CN-cycle are written as 



ET.CN 



^i+e^ + ei + ee + ^aiSi^ ^ ecN, (A5) 

(A6) 

The reciprocal of the time scale for collisional destruction 
of species i is 

K, = Ny^C,. (A7) 
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Fig. 8. Upper: Reciprocal of tlie timescale for collisional de- 
struction of ^^C, ^'^C, ^*N, and ^^N. Lower: Temperature de- 
pendence of nuclear energy generation by the CN-cycle. For 
a comparison, the result by Kawaler (1988) is shown with the 
dotted line. 

Since the C^'s depend only on temperature, and N can 
be expressed as a function of the mass density, p, the 
Lagrangian variations of A'i's are expressed as 

bKi bp , 5T 



P T 



(A8) 



where = (91nCi/91nT)p,j, is the temperature depen- 
dence of Ci . Taking the Lagrangian variation of equations 
(A1)-(A3), we have 



ia + 
ia + Xi 



SP 



(a._i + l)^ + (/3,_i+i^,;_i) 

P T 



(z = l,4) (A9) 
5T^ 



(i = 2,5) (AlO) 



[ai-i-l) V{Pz-i-Vi) — 

P T , 



(i = 3,6) (All) 
Then, we obtain the recursive formulae for the cti's: 



as 



ia + Ki 
K3 



ae, a.2 



ia + A; 



■(ai + 1). 



ia + K 



— (0:2-1), 0^4 



ia -I- 



as, 



A. 



"5 



ia + A5 
and for the /^^'s: 



(q!4 + 1), ae 



ia-^Ki 



■(«5-1) 



(A12) 
(A13) 
(A14) 



/3i 



ia -|- K\ 

ia + K:. 



(/36 + l'6-'^l), 132 



Xo 



ia + X2 



ia + K, 



(A15) 



ia + A5 ia + Kq 

(A17) 

Evaluating the at's and /3i's by equations (A12)-(A17), 
we obtain £p_cN and Et^cn from equations (A5) and (A6), 
respectively. 
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